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Abstract 

The theory of growth kinetics developed previously is extended to the asym- 
metric case of off-critical quenches for systems with a conserved scalar order 
parameter. In this instance the new parameter M, the average global value of 



the order parameter, enters the theory. For M = one has critical quenches, 
q [ while for sufficiently large M one approaches the coexistence curve. For all 



M the theory supports a scaling solution for the order parameter correlation 



function with the Lifshitz-Slyozov- Wagner growth law L ~ t ' s . The theoret- 
ically determined scaling function depends only on the spatial dimensionality 
d and the parameter M, and is determined explicitly here in two and three 
dimensions. Near the coexistence curve oscillations in the scaling function 
are suppressed. The structure factor displays Porod's law Q~^ d+1 > behaviour 
at large scaled wavenumbers Q, and Q 4 behaviour at small scaled wavenum- 
bers, for all M. The peak in the structure factor widens as M increases and 
develops a significant tail for quenches near the coexistence curve. This is in 



qualitative agreement with simulations. 
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I. INTRODUCTION 

In previous work [[]] a theory based on a field-theoretic Langevin model was developed 
to treat the growth kinetics of a system with a conserved scalar order parameter for the 
case of symmetric or critical quenches. In this paper the lowest order version of this theory 
is extended to off-critical quenches. Quenches to a final state near the coexistence curve 
where the volume fraction of the minority phase is small have been studied by a variety of 
approaches, but theoretical studies of the Langevin model have never been extended into 
this regime. The techniques developed thus far are generalizations of the Lifshitz-Slyosov- 
Wagner (LSW) treatment |2||| which considers one spherical droplet interacting with a mean 
concentration field. This approach is valid only in the limit of zero volume fraction, but other 
mean field theoretic and statistical mechanical techniques have been developed to incorporate 
the effects of the interaction of other droplets and extend the theory to slightly larger volume 
fractions PHT2J . Another approach to the problem is to use numerical simulations in concert 
with a theory describing the concentration field around spherical droplets (essentially an 
electrostatics problem with moving boundary conditions) fl3| . Direct simulations of the 



Langevin equation exist for the off-critical case in two dimensions [|14||T5|| , but we are not 
aware of any such simulations for three dimensions [ |T6|j . 

The theory developed in |]| shows how one can solve some of the thorny problems asso- 
ciated with growth kinetics for the conserved order parameter (COP) case. The theory can 
be evaluated as a well-defined sequence of approximations with qualitative and quantitative 



improvement as one moves along this sequence |17fl . In this paper we limit ourselves to the 
lowest order approximation in this formalism. From the work in [l[] we know that there 
are substantial limitations associated with this approximation and these are discussed in 
some detail in section VI. However, it is also known from [|IJ that this approximation gives 
good results for the scaling function for correlations of the order parameter. We therefore 
concentrate on this quantity here. 

The new element in this work compared to p| is that the average value of the scalar 



order parameter ip is no longer zero: 

(^(R,t)) = M. (1.1) 

M is independent of R and t because of the statistical homogeneity of the system and the 
conservation law, respectively. 

The main results of this paper are that, as in the critical case, the theory supports a 
long-time scaling solution for the order parameter correlation function 

C(K,t) = (5ij)(R,t)8ip(0,t)) 

=* F K) (h2) 

where 5ip = ip — M, ipo is the magnitude of the ordered field in equilibrium, and L(t) is the 
characteristic length in the theory, the average size of the domains. For later convenience we 
define the normalized quantity M = M/ip , which ranges between —1 and 1. When M — ±1 
the system is at the coexistence curve. It is found that, for all M and for long times t after 
the quench, the growth law L ~ i 1 / 3 holds. For small scaled distances x (= |R|/L(£)) one is 
able to find a scaling solution of the form 

F(x) = 1 - M 2 - e- y2/2 ax(l + fox + ■ ■ ■) (1.3) 

where the parameter y is related to M by 

M = erf(y/V2). (1.4) 



Unlike the non-conserved order parameter (NCOP) case treated earlier JT8|, the coefficient 
fa is not zero and must be determined, along with a, as part of a non-linear eigenvalue 
problem, fo is found to be negative for M = and monotonically decreases as M — *■ 1. 
Thus, in this theory, the Tomita sum rule (fa = 0) [|19|] is strongly broken as one approaches 
the coexistence curve. This appears to be an important limitation of the current theory. 

The scaling function, F, has been determined explicitly in two and three dimensions by 
solving the non-linear eigenvalue problem mentioned above. The dependence of F on M is 



weak for small M. As M increases the first zero of F moves to larger scaled distances and 
the first minimum of F becomes shallower. As M increases further the oscillatory behaviour 
is suppressed and the predominant behaviour is that of decay, as predicted by our large y 
asymptotic analysis. For large M there are oscillations at large x whose wavelength increases 
as one approaches the coexistence curve. These oscillations preserve the conservation law 

fd d xF(x)=0 (1.5) 

despite the existence of the strong decay. 

The structure factor is the Fourier transform of the correlation function and one has 

C(q,i) = iJ 2 L d (t)F(Q) (1.6) 

where Q = qL is a scaled wavenumber. F(Q) is characterized by five parameters. Since 
F(Q) is a peaked quantity, the peak position Q ma x, and height F(Q max ), are of interest as 
functions of M. The full-width at half- maximum, measured in units of Q maX i is also relevant. 
The linear term in (|1.3|) leads to Porod's law [2(| for the large Q tail of the structure factor 



- ( Y (d+1) 

F(Q) = F(Q max )A P (M) -¥- , (1.7) 

\ ^iraax I 



while for small Q, F(Q) behaves as GT 



F(Q)=F(Q max )A 4 (M)(-^-\ . (1.8) 

\ W max ) 

This behaviour at small Q is a result of a conserved diffusive field existing away from the 
interfaces which mediates the interaction among the interfaces. Our analysis shows that 
both Qmax and F(Q max ) decrease to zero as M — > 1. The width of the peak increases 
slightly for small M, but then develops a significant tail as M — > 1 for intermediate values 
of Q near the base of the peak. Ap(M) is a decreasing function of M, approaching zero in 
a cusp asM-* 1. The coefficient A±(M) increases with increasing M, growing rapidly near 
M — 1. Damped oscillations are also seen in the structure factor around the Q~ A behaviour 
at intermediate values of Q, a result also seen by other investigators |22|,[l3| . 



When we compare our results with those of other workers we find good qualitative agree- 
ment for both F(x) and F(Q) as functions of M. There are some quantitative differences, 
though. We believe that the difference in the form of F(Q) is due to our low estimate for 
the coefficient A±, leading to a peak in F(Q) which is too narrow. It seems likely that the 
lack of quantitative agreement is associated with the breaking of the Tomita sum rule in the 
theory. On the other hand, this is the only theory which has led to a determination of A&, 
and one can hope that using higher order approximations will give improved results. 

In the next section the theory forming the basis for this work is outlined. In section 
III the averages are performed which are relevant to the off-critical case. The end result of 
these manipulations is a non-linear equation for the scaling function. Section IV looks at 
the various limiting cases of the theory: the small x, large x, small Q, and large y behaviour. 
Section V presents the numerical study of the equation for F(x) in two and three dimensions. 
Comparison of the results of this paper with results from other investigators is made in 
section VI. The paper concludes with some comments about future areas of research and 
improvements to the theory. 

II. THE MODEL 

The dynamics are modelled using a noiseless time-dependent Ginzburg-Landau equation 
for a conserved scalar order parameter ip with a non-zero average M: 

^ = DoVliV'm)) -V 2 M1)}- (2-1) 

V(ip) is a double- well potential with degenerate minima at ±Vo> Du ^ is otherwise unspecified 
since, as we shall see, our results are independent of the precise form of V |p3fl . Here, the 
convienient notation in which 1 represents (Ri, ti) is used. The D Vl factor in ( |2.1| ) ensures 
that the system has a conserved order parameter. Thermal sources of noise are neglected 
because it is assumed here that the quench is to zero temperature. Randomess enters the 
problem through the initial conditions where we assume that the initial values of Sip = ip — M 
are governed by a Gaussian probability distribution characterized by 
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(V>(R,0)) = M (2.2) 

(8ip(R, 0)<5^(R, 0)) = e 5(R - R'). (2.3) 

Our final results are independent of the amplitude e appearing in the initial distribution. 

A method for extracting the correlation functions from ( |2.1| ) is described in M. Here we 
will merely outline the salient points. The order parameter is written as 

V>(1) = cr[m{l)] - u(l) (2.4) 

where a is the equilibrium interfacial profile and u represents fluctuations about this ordered 
value, m is assumed to be a random field whose zeros correspond to the zeros of a, that 
is, to the positions of the interfaces. The nature and interpretation of m will be discussed 
below. In the NCOP case the fluctuating field u can be safely ignored but in the COP case 
it is the field u which couples distant interfaces by permitting currents of minority phase 
atoms to flow through the matrix. In |l| the theory was closed by relating u back to a and 
m via the equation 

u(l) = ^<r(l) + AV?m(l) + --- (2-5) 

where Uq and A are parameters. This form satisfies the desired properties that u is conserved 
in bulk, odd in m, and 0(1/ L) everywhere. This last requirement comes from the fact that 
the interfaces are a source of u with a contribution proportional to the local curvature of an 
interface k (u\s ~ k). It can then be shown that if (|2.5| ) holds, u satisfies the familiar form 
V 2 m = away from the interfaces, a is chosen to satisfy the equation for an equilibrium 
interface 

\a 2 (l) = V'(a(m(l))) (2.6) 

where the factor 1/2 is inserted for convenience and cr n (m) = d n a(m)/dm n . The bound- 
ary condition, lim m ^ ±00 a = ±^o, guarantees that the system orders at the appropriate 
equilibrium value of the order parameter and results in the useful relation 



dxa x {x) = 2ip . (2.7) 



As shown in |1[, equations (|2.5| ) and ( p.6|) can be substituted into ( |2.1|) and the result 



multiplied by <r(2) and averaged to get an equation for 

C aa (R,t) = (a(l)a(2)). (2.8) 

At late times 1/L is expected to be small and one finds to leading order in 1/L: 

~C w (R,t) = -D q 2 V 2 (^-C aa (R,t) + AV 2 C CTm (R,t)) (2.9) 

where, q* = (V"(a)) = V"(if> ) + 0(1/ L) and C am (R,£) = (a(l)m(2)). Here, equal times 
are considered and statistical homogeneity of the system has been assumed so t\ = t 2 = t 
and R = R2 — Ri- Since 

C(R,t) = (#(1)^(2)) 

= <or(l)flr(2)> - M 2 + O Q) (2.10) 

we see that ( |2.9| ) is an equation for C(R, i) to 0(1/ L) since the action of the derivatives 
eliminates the disconnected part of the correlation function. 

A key aspect of the theory is the choice of the probability distribution P[m] governing 
the field m. This point is discussed in some detail in references |j],17|24 . Here we limit 



ourselves to the case where P[m] is given by an off-set Gaussian with 

m(t) = (m(l)) , (2.11) 

and (. . .)o is over a probability distribution Po[^ m ] which is Gaussian with respect to 5m(l) = 
m(l) — m(t). P()[5m\ is then determined by the variance 

C (12) = (Sm(l)8m(2)) . (2.12) 

Since the field m can be approximately interpreted as the perpendicular distance to the 
nearest interface, the off-set corresponds to a greater probability to be in one phase than in 
the other. The effects of this non-zero average will be explored in the next section. 
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III. EVALUATION OF AVERAGES: THE SCALING EQUATION OF MOTION 

A. Evaluation of Averages 



In order for ( |2.9| ) to be a closed equation for C^, it is necessary to relate C am to C a(7 . 



As in [0], this is done by using the Gaussian nature of m. Now, however, the Gaussian 
probability distribution must satisfy the condition (m(l)) = rh(t) ^ 0. Taking this into 
account one finds, using the standard properties of Gaussian integrals [0 , 



C CTm (12) = (or 1 (l)) C (12) + Mm{t). (3.1) 

The spatially independent term is eliminated by the action of the Laplacian in ( f2.9[ ). Since 
m is a Gaussian random field it follows that 

(a 1 (l)) = Jdx l( r 1 (x 1 )^x 1 ) (3.2) 

where 

$( Xl ) = l e -\{x,-rh(t)?/S (t) ,3 3) 

y/2nS (t) 
with 

S (t) = Co(ll). (3.4) 

Since m is a measure of the distance away from an interface it is expected that in the 
long-time scaling regime Sq ~ L 2 and m(t) ~ L. Therefore the limit 

m(t) 
lim . V ; = y (3.5) 



exists. In evaluating ( pT2| ) it is important to note that, for a wide class of potentials, 0i(xi) 
goes exponentially to zero for large \x\\. Therefore one can, after eliminating rh in favour of 



yJSo(t), expand ( |3.2| ) in powers of 5*0 1 and use ( p.7| ) to obtain 

2^ 



(a 1 (l)) = r " e-^ + O(S^). (3.6) 

/27T5 (t) 



The parameter y can be related to M by using the derivative relation |]25 

dM d(a) 



dm(t) dm(t) 
from which follows 



(<ri)o, (3-7) 



rm(t) / 2 _i * 2 , s 

M= / dzJ—j-r^e 2 ^>. (3.8) 

Thus, with M = M/ip , 

M = erf (y/V2) . (3.9) 

Therefore, there is a one-to-one correspondence between y and M. y = is a critical quench 
and y — > oo corresponds to a quench at the coexistence curve. 
With the definition 

/(R.«) = ^ (3-10) 

and the use of (|3.1|) and (|3.6|), the equation of motion (|2.9|) takes the form 



IJ^(R,t) = -Ax? 2 V 2 I ^C w (R,t) + xJ^ e-^ 2 V 2 f(R,t) ] . (3.11) 



The theory can be closed by relating Cq(R, t) to C CTCT (R, t) via the relation [18 



where C nm (R, t) = (<7 n (l)<7 m (2)) . The fact that m is Gaussian enables one to write 

Cn(R, £)= / dxidx2<7i{x\)ai{x2)Q{x\,X2) (3.13) 

where 

$( x x ) = 7 e ~ 5 ^ [( g i ~^W) 2 + (^2 -m(t)) 2 -2/Qri -m(t))(x 2 -m(t))] /g ^ 

27rS , 



with 7 = (1 — / 2 ) 2. Again, m(t) can be eliminated in favour of yJSo(t) and the result 
( P-13| ) expanded in terms of Sq 1 . At long times the leading order term dominates, resulting 



in 
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Cn(R,t) = ^re^7 f dx 1 dx 2 a 1 (x 1 )a 1 (x 2 ) + O(S 3/2 ) (3.15) 



where the integral is easily found to give 4?/>q. Integration of fl3.12|) with the substitution 



Q3.15J ), keeping in mind the definition of /, gives the desired relation between C aa and /: 

C„.M/'^L_ eT £. (3 . 16 ) 

VT JO Vl — S 2 

For a critical quench y — 0, and this reduces to the now standard expression C aa = 
fao sin -1 /. Equations ( |3.16| ) and (|3.11 ) form a closed system of equations for C aa . 



B. The Scaling Regime 

At long times the correlation function is expected to have a scaling solution. The ansatz 

a CT (R,t)=^ 2 F(x), (3.17) 



with x = |R|/L(i) being the scaled distance, can be substituted into equations ( [3.1 1|) and 
( ETI6I ) to give 

xF' = V 2 (uF + V 2 f) (3.18) 

2 ff ds i„u-i 

F = - ^= e -,y— (3.19) 

vr h Vl - s 2 

where 

F(x) = e y2/2 F{x) (3.20) 

u = 2UoDoq ° (3.21) 



and we have chosen L = [\u/iIjqUq)J2Sq/t^. Note that, for scaling we require that u is a 
constant and therefore L ~ t 1 ' 3 . One expects F to have a weaker y dependence than F for 
large y. It is useful to express / in terms of 



* = VlTT < 3 - 22 > 



in equation (|3.19|) to obtain the well-behaved integral representation 
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f = zLt-t^- ( 3 - 23 ) 



For future reference, the spatial derivative of ( |3.23| ) is 

dF _ 4 e" y2 * 2 / 2 <9$ 
dx n 1 + $ 2 dx 

It is noteworthy that F has a lower bound. Let $ — > oo in ( |3.23| ) and write 



(3.24) 



-^L^-^ A =y^ <^> 



I(y) has the following properties: 



w - -\ 



d n 

Ty I(y) = 2 

Integrating these equations yields 



1 — erf —p= 
\V2, 



4-er/(-?=). (3.26) 

dy JK V2 v ; 



2 



I(y) = -^[l-erfl [ ^=jj =-^(l-|M|) 2 . (3.27) 

Thus, the lower bound on F is 

F mm = -e y2/2 {l-\M\f. (3.28) 

AsM->l this lower bound approaches zero. 

Equations ( 3.18Q and ( p.23|) constitute a non-linear eigenvalue problem for F, with a 



unique solution determined by the boundary conditions at small x and the physical condition 
F — > exponentially as x — > oo. We will see that u is determined as part of the solution. 
The only parameters entering into the determination of F are y and the dimension d, which 
appears in the spherically symmetric Laplacian. 

IV. LIMITING CASES 
A. Small x behaviour 

The small x behaviour of F can be determined analytically. We find that F has the form 
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F = F(0) - ax(l + (3 2 x 2 + fox 6 + ■ 



(4-1) 



Expanding $ for small x 



$ = $iX + § 2 x + 



(4.2) 



and using (|3.24|) to connect the power series for F and $ gives 



$1 



$5 



7TOJ 
T 

7ra/3 2 



(4.3) 
(4.4) 



Substitution of these results into fl3.18|) give, at 0(l/x) 

-3n 2 af3 2 (d + 1) 



u 



F(0) can be determined from (|3.23|) by noting that $(0) 
technique that was used to derive ( |3.2S| ) . The result is 

F(0) = (1 - M 2 )e y2/2 . 



(4.5) 
and then using the same 

(4.6) 



The equation of motion fl3.18| ) can be partially integrated using the Green's function for the 
Laplacian p5" |. For d > 2 the result is 



uF + V 2 / 



d-1 



d 



d-2 



X 



I d {x)-2(I 2 (x)-I 2 {oo)) 



(4.7) 



with 



I d (x) = / z d - 1 F(z)dz. 
Jo 



(4.8) 



The results ( POT ) and ( Of ) can be substituted into ( J4.7| ). Since one can show that 

/d(x) 



lim 



T d-2 x _»o 



a;->U x 



lim J 2 (x) = 



(4.9) 



one has at 0(1) 



2 r , s -K 2 a 2 d -, . 
I a (oo) = _+flF(0). 



d-2 



(4.10) 
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For d = 2 the analogous results are 

uF + V 2 / = (1 - 2 lnx)/ 2 (x) + 2(J 2 (x) - J 2 (oo)) (4.11) 

with 

J 2 (x) = r^ zF(z)ln;z. (4.12) 



/o 
Since lim^o ^2 (a?) = substitution of ( |4.3| ) and (|4.4|) into ( [4.1 1| ) give, at C(l) 



7T 2 « 2 



-2J 2 (oo) = — + uF. (4.13) 

Knowledge of the parameters a and /3 2 allows one to determine the constants u and I 2 (oo) 
(or J 2 (oo)) appearing in the equation of motion for a given value of y. Numerically, what 
this means is that values for a and /3 2 are chosen and equations (|4.7| ) (or (|4.11|) ) and Q3.23|) 



integrated forward in x. a and /3 2 are adjusted until F satisfies both the conservation law 
([L.5D and the boundary condition F — > exponentially as x — >• 00. 

B. Large x behaviour 

For large x both F and / are small and $ can be expanded about its asymptotic value 

§(x) = l + ri(x). (4.14) 

Substitution into (|3.24| ) and integration yields a relation between F and rj 



F(x) = --e- y2/2 r](x). (4.15) 

7T 

Also, from ( |3.22| ) we have f(x) = —r](x) so we may rewrite ( |3.18| ) in the large-x limit as 

x p = V 2 {uF + -e y2/2 V 2 F). (4.16) 

Li 

This equation supports damped oscillatory solutions of the form 

F ~ F x--/»exp (-^7;^) cos (^^ + ^ (4 . 17) 

where Fq and are constants. Note that as y — >• 00 the wavelength of the oscillations 
increases and the exponential term goes to 1. This means that in this limit one must go to 
progressively larger values of x before one sees this asymptotic behaviour. 
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C. Small Q behaviour 

It is the Fourier transform of the order parameter correlation function, 

C(q,t) = f d d Re i * n {i/j(R,t)'ip(0,t)) 

= 4> 2 L d (t)F(Q) + (2n) d M 2 5(q) (4.18) 

with Q = qL(t) and it is 

F(Q) = J d d xe iQ *F(x), (4.19) 

that is measured in a scattering experiment. In the total scattering cross-section we expect 
that at long times there is a dynamic contribution to the forward Bragg peak, 

)imiftL d {t)F(Q)=A5{q), (4.20) 

in addition to the static contribution (2ir) d M 2 5(q). Since 

A= fd d q \im *p 2 L d {t)F{Q) 

J t— >oo 

= (2tt)> 2 F(0) = (270^(1 - M 2 ), (4.21) 

the total contribution to the forward Bragg peak at late times is (27r) d ipQ5(q), as expected. 
To examine the small Q behaviour of the structure factor F(Q), it is useful to consider 
the moments of F(x) 

W p = fd d xx p F(x) (4.22) 

which can be found by multiplying (|3.18|) by x p and integrating. The result is that Wo = 
W 2 = while 

U{d + 2)^ y2 f. 



Thus we have 



W 4 = y ' e~ y / d d xf(x). (4.23) 



F(Q) ~ AQ A (4.24) 
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for small Q where 



A = - JTI ld d xf(x). (4.25) 



The Q 4 behaviour of F(Q) at small Q is a consequence of the fact that in the theory u 
is conserved away from the interfaces, and this behaviour does not depend on the specific 
ansatz for u. In order for F(Q) to be positive definite ( |4.25| ) implies that 

d d xf(x) < 0. (4.26) 



This means that lim q ^ {\m q (t)\ ) / S < 0. As pointed out by Yeung et al. |26| this is a 
shortcoming of the fact that m is a Gaussian variable. This problem is resolved when 
non- Gaussian corrections to the probabilty distribution for the field m are considered [|T|. 

D. Large y behaviour 

An analytic result for the limit as one approaches the coexistence curve, M — > 1 , y — > oo 
is of interest because it allows one to make comparisions with other theories developed for 
this regime. From the numerical analysis in the next section the following facts emerge. The 
first is that, as y increases, the scaled length x over which the correlation function takes 
significant values decreases. This suggests that x should be rescaled as 

x = y p z (4.27) 

with p < 0. Second, it appears that u grows as some power of y for large y. We are led to 
assume the form 

u = u^y" (4.28) 

with n > 1. Finally, since the interesting behaviour occurs near the origin where the quantity 
$ is small, we can rescale $ 



M 



$ = - < 1 (4.29) 

y 
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where grows slowly but does not break the bound near the origin. Using this definition in 
equation ( |3.22| ) leads to 



/«l-2- 



(4.30) 



Armed with these results we proceed to re-examine the theory. Using (|4.29 ) in equation 
( |3.23| ) and letting t = ys in the integrand gives 

dt 



4 rv 
F{z) = - 



-t 2 /2 



Try Jftz) 1 + (t/y)' 



To leading order in y 1 this is 



where 



F(z) = — / dte-' 2 ' 2 = -F^z) 
ivy Jftz) y 



'-'-■'-''^-^(vf 



Under the rescaling outlined above the equation of motion (|4.7|) for d > 2 becomes 



Uooy^F^z) - 2y 



-2(p+l)v72 ,2 _ V 



2p-l 



d-2 



yd-2 



I d (z)-2{I 2 (z)-I 2 (oo)) 



with 



h(z) 



dss^F^is). 



(4.31) 



(4.32) 



(4.33) 



(4.34) 



(4.35) 



Since p < the right hand side does not contribute as y — » oo. The integrals in the d = 2 
case also do not contribute in this limit so the following results are valid for d > 2. Balancing 
powers of y on the left hand side of Q4.34 ) gives 

n+1 



V 



and 



Ury^r r-yQ Z V 



2 J. 2 



(4.36) 



(4.37) 



Using ( |4.33| ) this becomes a simple equation for 



"^V- erf \T2, 



72/2 



V z d) z . (4.38) 

In the numerical solution of ( |4.38| ), Uoo is a parameter which is found from a fit of the 



numerically determined u to the form (|4.28j ) for large y. 
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E. Small and large z behaviour for large y 



An examination of ( [4.38 ) in the limit of large and small z is instructive. When z is small 



is expected to be small so can be expanded as a power series in z 



z + (j) 2 z 2 + ■■■. (4.39) 



Also, in this limit, ( |4.38| ) simplifies to the form 

*2 ( [2 



Uoo\l- [ 1- W-0 | =VV. (4.40) 



Matching powers of z to leading order gives 



1 \ix) \2d 



(4.41) 



This result allows one to make predictions about the asymptotic behaviour of a and fa when 
y is large. Using (|4.27| ) and Q4.33|) one can write 

F(x) = -\---(f) 1 y- p - 1 x + --- (4.42) 

y V 7T n 

giving 

a = a^- 1 (4.43) 

for large y with a^ = Afa/ir. If we also write 

fa = P 2 ooy m (4.44) 

then 

_ = _3AA*l) r ,, = _^„ (4 4g) 



Matching the coefficient and the exponent leads, using fl4.36| ), to the following relationships: 

a = «ooi/ (n - 1)/2 (4.46) 

& = /W n+1)/2 (4-47) 

«oo = -- ^= f ^ &OC. (4.48) 

V27r \ a 



Thus, a graph of fa /a vs. y for large y will be linear with a slope that depends only on the 
dimensionality of the system. 

At large z, <j) is large and ( [4.38 ) is well approximated by 

u^—e-^l 2 = VV- (4.49) 

ir<p 

For d > 2, standard asymptotic analysis yields, at next to leading order: 
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This implies that 



:;, = taVa lsi(sfe)l- < 450 » 



- , x 4(d-2) 1 . 

Foo{z) = -A_ >-- (4.51) 

Uqq Z 

for large z and d > 2. For d = 2 one has 

FcJz) = ^ . (4.52) 

It is clear that, near the coexistence curve, the oscillations in the scaling form become 
insignificant and are dominated by a strong decay. 

V. NUMERICAL SOLUTION OF THE NON-LINEAR EIGENVALUE PROBLEM 

A. a and fa as a function of y 



In this section the numerical solution of (|4.7|) (or (|4.11|) ) coupled with ( |3.23|) in two 



and three dimensions will be discussed. The equations are integrated forward from x = 
using a fourth order Runge-Kutta integrator with step size Sx = 0.001, subject to the initial 
conditions F(0) = (1-M 2 )e^ 2 , F'(0) = -a , F"(0) = -2ap 2 , and f (0) = 0. This method 
of integration seems numerically stable and insensitive to the choice of Sx. The search for the 
eigenvalues a and fa involves requiring that the solution F obey the conservation law ( |1.5| ) 
and have the physically acceptable behaviour F — > exponentially as x — > oo. This search 
is performed by fixing a and then searching for the value of fa which pushes the diverging, 
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unphysical solution to larger values of x. The value of a is then adjusted so that the flat 
region of F at large x is properly zeroed. The procedure is repeated with the new value of 
a until the exponentially growing solution is pushed as far from the origin as possible and 
until F is zeroed as well as possible. The degree to which the conservation law is satisfied 
naturally depends on how well the function is zeroed. The convergence of the eigenvalues is 
fast and F can be zeroed to better than 10 -6 using this method. 

The results for the eigenvalues a and fa are shown in Fig. |l|. One sees that a initially 
decreases reaching a minimum at y ~ 1 and then rapidly becomes large and positive as the 
coexistence curve is approached. The eigenvalue fa is negative at y = and monotonically 
decreases as y increases, decreasing rapidly as y becomes large (M — >■ 1). Equation (|4.48|) 
predicts that a graph of fa/a will be linear for large y with a slope -0.157 for d = 3 and 
-0.139 for d = 2 . We see this linear behaviour and have measured slopes of about -0.143 and 
-0.129 for d = 2 and 3 respectively. The exponents n and p, and the coefficient Mqo defined 
in (|4.27|) and (|4.28 ) can be found by fitting the large y behaviour of a and fa to the forms 



( |4.46| ) and ( (4.47[ ). For three dimensions one finds n ~ 8, p ~ —4.5 and Mqo ~ 0.0033. In 
two dimensions one has n ~ 6, p ~ —3.5 and u^ ~ 0.024. When considering these results it 
should be kept in mind that only a few values of y around y = 4 were used to obtain these 
values. In principle, both the exponents and the coefficient can be accurately obtained by 
extending the numerical analysis to larger values of y. In practice, this is difficult due to 
reasons that are discussed below. 

B. Scaling Function as a function of M 

The dependence of the scaling function F(x) on M is shown in Fig. |2| for two dimensions 
and in Fig. |3] for three dimensions. In these plots F(x) is normalized so that F(0) = 1. 
Both figures show that F(x) depends only weakly on M for values of M < 0.4. The scaling 
function has a prominent oscillatory component which is necessary to satisfy the conservation 
law. At intermediate values of M, the position of the first minimum of F(x) occurs at larger 
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values of x and the depth of this minimum decreases as M increases. The depth of the 
oscillations is greater in two dimensions than in three. These stronger oscillations make the 
presence of the lower bound on F(x) noticeable, and near the coexistence curve the minima 
in the scaling function are very flat in order to be consistent with this bound. 

As M — > 1 the scaling function approaches its asymptotic form ( |4.33| ) which can be 
determined by numerically solving (|4.38|) using the values of u^ found in the previous 
section. Since we know the exponent p, we can rescale the distance x using ( |4.27|) and plot 



Foo(z) = yF(xy p ) for large values of y. This is done for two and three dimensions in Fig. |j. 
The asymptotic forms obtained by solving (|4.38|) are also shown in this figure. We see that 



F(x) decays very rapidly when the system is near the coexistence curve. Oscillations do 
occur for these values of y, but they occur at large x and have a small amplitude and large 
wavelength. For d = 3 the curves appear to approach the asymptotic form as y increases. 
Hovever, for d = 2 the asymptotic form is not approached if one uses Mqo = 0.024. A value of 
Uoo = 0.036 gives a better fit and the form obtained using this value is the one shown in Fig. 
[|. Matching the asymptotic form to the rescaled large y scaling function is another way to 
determine Uoo. We believe that the two methods for finding Uoo give different values because 
in the fit of a and fa to the forms (|4.46|) and ( |4.47| ) we do not have values of y which are 



large enough to be in the asymptotic regime. Larger values of y are difficult to reach because 
one runs into numerical problems as the theoretical lower bound on F(x) approaches zero. 
These numerical problems are especially significant in two dimensions since the oscillations 
in the correlation function are stronger than in three dimensions. 

C. Scaling of the Structure Factor 

The structure factor, 

F(Q) = f d d xe iQ *F{x) (5.1) 

was calculated by taking the Fourier transform of our numerically determined F(x). We 
find that as M increases the height of the peak decreases and the peak position moves to 
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smaller values of Q. Graphs of the normalized structure factor for various M are shown in 
Fig. |5] for d = 2 and Fig. |6] for d = 3. Logarithmic plots reveal the power- law dependence 
of F(Q) for large and small Q (Fig. |7|). For small Q, F(Q) ~ Q 4 in both d = 2 and 3, for 
all M. Small deviations from the Q A behaviour can be seen, but we atttribute these to the 
unreliablity of the numerical determination of F(x) for extremely large x. For large Q, one 
observes Porod's law, F(Q) ~ Q~~( d+1 \ for all M. The coefficients A^ and Ap defined in the 
introduction are determined and plotted in Figs. |8] and § respectively as functions of M. A4 
increases with increasing M and Ap is a decreasing function of M, approaching zero like a 
cusp at the coexistence curve. 

Figures [5] and |6] show that the width of the peak increases slightly, but is rather insensitive 
to changes in M until very near the coexistence curve. In the logarithmic plots there 
appear to be damped oscillations in F(Q) at intermediate Q before the onset of the Q^^ 1 ^ 
behaviour. In two dimensions, as one approaches the coexistence curve, the main peak 
decreases in amplitude until it is comparable to these oscillations, which show up as a 
shoulder to the main peak. In three dimensions there is a tail on the large Q side of the 
peak in the structure factor, which also grows as M increases. Both the secondary peak and 
the tail may be related to the fact that the Tomita sum rule is stongly broken as M — > 1. 
The large coefficent of x 2 in the small x expansion of F(x) will lead to corrections to Porod's 
law for the medium Q behaviour of the structure factor. 

VI. COMPARISONS 

In order to test the validity of the assumptions made in this paper the results for F(x) 
and F(Q) will be compared with the relevant results of other investigators. Experiments 
involving neutron scattering off of a binary alloy have been done for a fixed M p7|| , but 
we have been unable to find any experimental study of the dependence of F(Q) on M. 
One problem with doing experiments near the coexistence curve is that the small volume 
fraction of the minority phase causes the structure factor to have a small amplitude, thus 



22 



making it difficult to measure. Another problem that arises when comparing experiment to 
theory is that it is unclear what volume fraction was used in a given experiment, making a 
straightforward comparision difficult. 

While there are high quality numerical simulations for critical quenches |28|^3~0f , there has 
been far less work on off-critical quenches. One example is the direct numerical simulation of 
the Cahn-Hilliard equation in two dimensions performed by Chakrabarti et al. [y|. Here, we 
compare their result for the correlation function with ours. Their functions are scaled so that 
the first zero of the correlation function occurs at x = 1, and we have adjusted our length 
scale to correspond to this. The comparisions for volume fractions <fi = 0.5, 0.21 and 0.05 
are shown in Fig. |I| The relationship between M and the volume fraction is 

0=I(1-M) (6.1) 

which is valid for quenches to T = 0. The quantitative agreement is poor. In particular, 
the theory predicts that at large x the oscillations in F(x) have much larger amplitude than 
seen in the simulations. Nevertheless, the positions of peaks and troughs of the oscillations 
are in qualitative agreement. In addition, we agree on the observation that oscillations in 
the correlation function become weaker and have longer wavelength as the coexistence curve 
is approached. In summary, the qualitative agreement is reasonable. We are unaware of any 
direct simulation of the Cahn-Hilliard equation in three dimensions for the off-critical case. 
Such simulations are difficult because large system sizes are required to give a statistically 
meaningful distribution of droplets when the volume fraction is small. 

One can also make comparisons with generalizations of the LSW theory ||4^[P2]j . This is 



in the regime of Ostwald ripening pi| . While much of the analysis in this case has focussed 
on the droplet distribution function, more recently a number of authors have determined 
F(Q) . In particular, here we will compare our three dimensional results with those of Akaiwa 
and Voorhees ||13|| . They assume that the droplets are spherical and interacting essentially 
electrostatically through a concentration field with both monopole and dipole contributions. 
Both the droplet size distribution function and the structure factor can be extracted by 
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numerical simulation of the equations produced by the theory. Our structure factors and 
those of [13|] are compared in Fig. [11]. Both results agree and give Porod's law at large 
Q. At small Q both results exhibit Q 4 behaviour, although our results seem to have a 
smaller coefficient of Q A than theirs. This may also be why the widths of our peaks are 
consistently smaller than those of [p|. There is also significant disagreement on the shape 



of the structure factor for values of Q just above the peak. In the theory presented here this 
regime of Q may be strongly affected by the breaking of the Tomita sum rule. 

VII. CONCLUSION 

In this paper it has been shown that the theory developed in |l[ can be extended to the 
case of off-critical quenches. The LSW t 1 ^ 3 law and the associated scaling behaviour are 
determined for the entire concentration range. The scaling function is a function only of 
the parameters d and M, changing significantly only close to the coexistence curve where 
the oscillations observed in the critical case are damped out. The structure factor exhibits 
Porod's law for large Q and Q A behaviour at small Q. This is the first theory which is 
capable of sensibly treating spinodal decomposition over the entire concentration range. 

As discussed above, there are a number of virtues of this theory. However, there are also 
important limitations. First, we have not been able to make contact with the LSW theory in 
the M — > 1 limit. This will require extending the current theory (or some improved version) 
to treat the droplet distribution in the dilute limit. This is a difficult but, to us, interesting 
challenge. Secondly, it is clear that we must extend the theory developed here to include non- 
Gaussian corrections if we are to remedy the problem of Co(q,t) going negative for small 
q. Since one expects this quantity to enter the determination of the droplet distribution 
function in an important way it is crucial to include non-Gaussian corrections if one is to 
make progress in this area. Non-Gaussian corrections have already been used to treat the 
critical COP case and this is discussed in ]I|. Finally, it seems reasonable to assume that 
the primary reason that we do not obtain good quantitative agreement for F{Q) and F(x) 
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is that we do not satisfy the Tomita sum rule. We speculate that in order to satisfy the 
Tomita sum rule an improved treatment of the gradient term in the consituitive relation 
75]) is required. 
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FIGURES 
FIG. 1. The eigenvalues a and /% as a function of y. (a) a for two dimensions (upper curve) 

and threedimensions (lower curve), (b) /% for two dimensions (lower curve at y = 0) and three 

dimensions (upper curve at y = 0). 

FIG. 2. The normalized scaling function F{x) in two dimensions for various M. In terms of 
decreasing depth of the first minimum the curves correspond to (a) M = 0, 0.2, and 0.4. (b) M = 
0.6, 0.8, and 0.9. 

FIG. 3. The normalized scaling function F{x) in three dimensions for various M. In terms of 
decreasing depth of the first minimum the curves correspond to (a) M = 0, 0.2, and 0.4. (b) M = 
0.6, 0.8, and 0.9. 

FIG. 4. The large y asymptotic scaling function F^z). From lowest to uppermost the solid 
curves correspond to y = 2.5, 3, 3.5, and 4 (M = 0.9876, 0.9973, .9995, and .9999 respectively). 
The infinite y form (M = 1), obtained from our asymptotic analysis is shown as a dashed line, (a) 
two dimensions, (b) three dimensions. 

FIG. 5. The normalized structure factor in two dimensions. From lowest to uppermost the 
curves correspond to M = 0, 0.4, 0.6, 0.8, and 0.9. 

FIG. 6. The normalized structure factor in three dimensions. From lowest to uppermost the 
curves correspond to M = 0, 0.4, 0.6, 0.8, and 0.9. 

FIG. 7. Logarithmic plots of e y > 2 F(Q). At ln(Q) = 0, from lowest to uppermost the curves 
correspond to M = 0, 0.4, 0.6, 0.8, and 0.9. (a) two dimensions, (b) three dimensions. In both 
graphs Q 4 behaviour is seen at small Q and Porod's law Q~( d+1 > behaviour occurs at large Q 
(dotted lines). 

FIG. 8. The coefficient A^{M) appearing in the small Q behaviour of F{Q) for two (upper 
curve) and three (lower curve) dimensions. 
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FIG. 9. The Porod's law coefficient Ap(M) for two (lower curve) and three (upper curve) 
dimensions. 

FIG. 10. A comparison of our scaling forms for the correlation function (solid lines) in two 
dimensions with those of |II| (dashed lines). The horizontal axis has been chosen so that the first 
zero of F{x) occurs at x = 1 for both functions, (a) (j) = 0.5. (b) <f> = 0.21. (c) <ft = 0.05. 

FIG. 11. A comparison of our structure factors (solid lines) in three dimensions with those of 
pi (dashed lines), (a) c/> = 0.3. (b) <j) = 0.2. (c) </> = 0.1. (d) <f> = 0.05. 
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